A lattice study of the Faddeev-Niemi effective action* 
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We perform a lattice analysis of the Faddeev-Niemi effective action conjectured to describe the low energy 
sector of SU(2) Yang-Mills theory. We generalize the effective action such that it contains all operators built from 
a unit color vector field n with 0(3) symmetry and maximally four derivatives. To avoid the presence of Goldstone 
bosons, we include explicit symmetry breaking terms parametrized by an external field h of mass-dimension two. 
We find a mass gap of the order of 1.5 GeV. 



1. Introduction 

Recently, Faddeev and Niemi (FN) have sug- 
gested that the infrared sector of Yang-Mills the- 
ory might be described by the following low- 
energy effective action M, 



S'fn = d x 



m 2 (5 At n) 



4e 2 



(n • <9 M n x d„n) 2 



(1) 



Here, n is a unit vector field with values on S , 
n 2 = n a n a = 1, a = 1,2,3; m 2 is a dimensionful 
and e a dimensionless coupling constant. The FN 
"field strength" is defined as 



(2) 



FN claim that (|l|) "is the unique local and 
Lorentz-invariant action for the unit vector n 
which is at most quadratic in time derivatives so 
that it admits a Hamiltonian interpretation and 
involves all such terms that are either relevant or 
marginal in the infrared limit" . 

It has been shown that Sfn supports string- 
like knot solitons characterized by a topo- 
logical charge which equals the Hopf index of the 
map n : S 3 — ► S 2 . In analogy with the Skyrme 
model, the H 2 term is needed for stabilization. 
The knot solitons can possibly be identified with 
gluonic flux tubes and are thus conjectured to cor- 
respond to glueballs. For a rewriting in terms of 
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curvature free SU(2) gauge fields and the corre- 
sponding reinterpretation of Sfn we refer to H . 

In this contribution we are going to address the 
following problems: First of all, neither the inter- 
pretation of n nor its relation to Yang-Mills the- 
ory have been clarified. An analytic derivation 
of the FN action requires an appropriate change 
of variables, A — > (n, X), which decomposes the 
Yang-Mills potential A into (a function of) n and 
some remainder X . Although progress in this di- 
rection has been made (| [j| , there are no conclu- 
sive results up to now. 

Second, there is no reason why both operators 
in the FN "Skyrme term" , which can be rewritten 
as 



H 2 = (9 M n) 4 + (0„n • cU) 2 



(3) 



should have the same coupling. Third, and con- 
ceptually most important, S'fn bas the same 
spontaneous symmetry breaking pattern as the 
non-linear cr-model, SU(2) — > £7(1). Hence, it 
should admit two Goldstone bosons and one ex- 
pects to find no mass gap. 

We have scrutinized the FN action using lat- 
tice methods. To this end we made a sufficiently 
general ansatz for an n-field action that contains 
(Q) as a special case. In particular, we allow for 
explicit symmetry breaking terms to avoid the ap- 
pearance of Goldstone bosons. 
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2. Method 



i 



After generating SU(2) lattice configurations 
using the standard Wilson action we fix to a co- 
variant gauge [fffi. We chose the Landau gauge 
(LG) defined by maximizing „ tr n L7^(x) 
w.r.t. the gauge transformation Q, leaving a 
residual global SU (2)-symmetry. The field n 
is then obtained via maximizing the functional 
Fmag = J2 X n tr ( T 3 9 U^ (x)r 3 9 Un(x)) of the max- 
imally Abelian gauge (MAG) fTjjfn]]. This yields 
a gauge transformation g which we use to define 
our n-field, 



(4) 



ft is important to note that this definition leaves 
a residual local U(l) unfixed. 

Since the configurations generated originally 
are randomly distributed along their orbits, the 
gauge fixing is absolutely crucial for rendering the 
definition almost gauge invariant Egj . 

Our ansatz for the effective action is S c g = 
Yli \Sj\n\ with couplings Xi and operators Sj. 
Up to fourth order in a gradient expansion there 
are the symmetric terms 



(c^n) 2 , (Dn) 2 , (9 M n) 4 , (d^n ■ 3„n) 2 , 



(5) 



and the symmetry breaking terms including a 
"source field" h, 



n • h , (n • h) 2 , (d„n) 2 n • h . 



(6) 



The couplings can be obtained by use of an 
inverse Monte Carlo method |l|] , where the (bro- 
ken) Ward identities for rotational symmetry pro- 
vide an overdetermined linear system, 
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(^ 6 [n]5f b [n,h])A J - = (7f[n]). 



(7) 



Here, F % ab and If are known functions of n, typi- 
cally linear combinations of n-point functions. 

All computations have been done on a 16 4 - 
lattice with Wilson coupling (3 — 2.35, lattice 
spacing 0.13 fm and periodic boundary condi- 
tions. For the LG we used Fourier accelerated 
steepest descent 14 . The MAG was achieved 



using two independent algorithms, one (Al) be- 
ing based on 'geometrical' iteration jy| , the other 
(All) analogous to LG fixing (see Fig. pi). 
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Figure 1. Behavior of the M AG-functional using 
different algorithms. 



3. Results 

As expected, we observe a non-vanishing ex- 
pectation value of the field (one-point function) 
in the three-direction that can be thought of as 
a 'magnetization' SDT, (n a ) = M6 a3 . Thus, the 
global symmetry is broken explicitely according 
to the pattern SU{2) — > U(l). This also shows 
up in the behavior of the two-point functions 
(Fig. 0), which exhibit clustering, (n 3 (0)n 3 (x)} ~ 
(n 3 )(n 3 ) = 9Jt 2 , for large distances. Furthermore, 
the transverse correlation function (of the would- 
be Goldstone bosons) 

G-^s) = ^(OyOr)}, i = l,2, (8) 

decays exponentially as shown in Fig. [| This 
means that there is a nonvanishing mass gap M 
whose value can be obtained by a fit to a cosh- 
function. 

The numerical values of the observables, 9Jt, 
M and the transverse susceptibility, x = 
^2 x G ± (x), are summarized in Table 1 for both 
algorithms: The slight disagreement between Al 
and All is expected from our still somewhat low 
statistics. The numerical results for the mass gap 
M lead to a value of about 1.5 GeV in physical 
units. The last column is a measure for the ac- 
curacy of the minimal ansatz consisting of the 
first (leading) terms of (||) and (||), respectively. 
In this case the (continuum) mass gap is deter- 
mined by the "source", M 2 = |h|. In addition, 
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Figure 2. Two-point correlators of the field n 
obtained via algorithm AI. The dotted line repre- 
sents the (squared) VEV of n, (n 3 ) 2 = 97t 2 . The 
same behavior is obtained via All with slightly 
different plateau value (see Table 1). 



Table I 

Numerical values for some observables (all num- 
bers in units of the lattice spacing) . 





Tt 


X ± 


M 


X^M 2 


AI 


0.436 


0.636 


0.95 


0.53 


All 


0.352 


0.596 


1.01 


0.58 



one has the exact Ward identity Wl = x^M 2 . Us- 
ing this relation one obtains the rough estimate 
that M ~ 1.2 GeV. Compared to the 'exact' (fit- 
ted) value of M ~ 1.5 GeV we find a qualitative 
agreement already to lowest order. 

The effective couplings have to be determined 
by solving (Q). Results already obtained will be 
reported elsewhere. 
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Figure 3. The transverse correlation function G , 
fitted to G ± (x) ~ cosh(-M(ai - L/2)). 
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